Standing gravitational waves from domain walls 
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We construct a plane symmetric, standing gravitational wave for a domain wall plus a massless 
scalar field. The scalar field can be associated with a fluid which has the properties of 'stiff' matter, 
i.e. matter in which the speed of sound equals the speed of light. Although domain walls are 
observationally ruled out in the present era the solution has interesting features which might shed 
light on the character of exact non-linear wave solutions to Einstein's equations. Additionally this 
solution may act as a template for higher dimensional 'brane-world' model standing waves. 

PACS numbers: 04.30.-w, 11.27.-|-d, 98.80.Cq 
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I. INTRODUCTION 

Most discussions of gravitational waves involve taking 
the weak field limit of Einstein's field equations thus ig- 
noring the non-linear character of the theory. There are 
examples of wave solutions to the full non-linear field 
equations such as the Einstein-Rosen cylindrical waves 
[l|, or non-linear gravitational plane wave solutions 0. 
These exact wave solutions to the full non-linear equa- 
tions are useful for understanding gravitational waves 
without the weak field assumption. 

In this article we present an exact standing gravita- 
tional wave solution in the presence of a domain wall 
and a massless scalar field. The solution we give here is 
a non-trivial modification of the usual domain wall solu- 
tion 



^Jl + k\z\ 



{dt^ - dz^)-{l+k\z\)[dx^ -^dy^] (1) 



i.e. it is not simply a weak field gravitational wave in the 
background ([T|), but is a solution to the full non- linear 
gravitational field equation. Unlike the pure domain wall 
solution, which is a vacuum solution everywhere except 
at the location of the wall i.e. z = 0, our solution is 
supported by a scalar field. 

Vacuum standing wave solutions to the Einstein field 
equations were investigated in Q and it was found that 
some special case of the Einstein- Rosen cylindrical waves 
may be considered standing waves. However because of 
the non-linearity of general relativity one can not (as 
in electrodynamics) take two oppositely traveling plane 
waves (i.e. the pp-wave solutions [3]) and add them to- 
gether to get a standing wave. Our plane symmetric. 
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standing wave solution is supported by a domain wall 
and scalar field. We find that the scalar field and the 
transverse metric components oscillate with a 7r/2 phase 
difference which suggests a passing of energy back and 
forth between the scalar and gravitational field. Since 
the energy of the scalar field can be localized this seems 
to imply that one can for this solution have a localized 
energy for the gravitational field - or at least for the wave 
part of the gravitational field. 

Defining a localized gravitational field energy is diffi- 
cult [5| , and strictly speaking it is impossible since via the 
equivalence principle one can always locally go to a frame 
in which space-time is Minkowski, i.e. a frame where the 
gravitational field vanishes. In the conclusion we will 
discuss in what limited sense it may be valid to think of 
localized gravitational field energy for this solution. 

The requirement of a domain wall in our solution is 
not surprising since even in electrodynamics one must 
have conducting planes in order to make standing waves. 
In the conclusion we argue that - minus the scalar field 
- there is a large similarity between the present solu- 
tion and the simple electro-magnetic standing wave be- 
tween two infinite conducting planes. For the solution in 
this paper there is only one plane - the domain wall at 
z = 0. The other 'wall' is provided by the static part 
of the gravitational field which increases as one moves 
away from z = 0. This is similar to the experiments 
of [6|] where neutrons were trapped between a refiecting 
surface at the surface of the Earth and the gravitational 
potential of the Earth. As in these neutron experiments 
we find discrete frequencies of oscillation for our gravi- 
tational standing wave just as the neutrons had discrete 
frequencies/energy levels. Also this analogy between the 
present gravitational solution and the neutrons exper- 
iments (where the oscillating part of the gravitational 
field is equivalent to the neutrons) again is suggestive 
that at least for the wave part of the gravitational field 
one has some idea of localization of gravitational energy. 
For neutrons there is a well defined concept of localiza- 
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tion of energy. 



II. METRIC AND FIELD EQUATIONS 

We begin by modifying the domain wall metric of ([T]) 
by adding t and z dependent ansatz functions S{t, z) and 
V{t,z) as follows 



Vl + kz 



-(1 + fcz) 



{dt^ - dz^) 



(2) 



Here fc > is a constant. Note that we write (1 + kz) 
instead of the standard domain wall function (l + fc|z|) so 
our solution applies only in the region z > 0. To obtain 
the equivalent solution for z < we could set fc < 0. If 
one sets k = and S — the form of the ansatz in ^ 
is similar to a special case of the traveling wave ansatz 
considered in [J Q . The ansatz functions of these lat- 
ter references are functions oiu — t — z or v — t + z as 
one would expect for traveling waves, whereas the ansatz 
functions of the metric in Q are simply functions of z 
and t without any special form, as one would expect for 
standing waves. With fc ^ and 5* 7^ one finds that 
the ansatz in ^ is some combination of the domain wall 
solution of ^ and the colliding plane wave solutions of 
p^ . The metric in ^ differs from the colliding plane 
wave solutions in that the domain wall function, {1 + kz), 
replaces the time-like coordinate, t, used in p^ . It is im- 
portant to stress the the standing wave solution we find 
below is not a linearized perturbative wave solution on 
some background, but is an exact solution. This solution 
is therefore mathematically distinct from situations such 
as, for example, the study of aspherical perturbed evo- 
lution of a relativistic star interacting with gravitational 
waves. 

In addition to the metric ansatz given above we take 
a massless scalar field, <j){t,z), which obeys the Klein- 
Gordon equation in the background from ^ 
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= VTTk~z (^4> - 0" - YTVz'^') " ° ' 

where overdots mean time derivatives and primes stand 
for derivatives with respect to z. (A system similar to 
that in ^ and ^ was recently investigated in terms of 
cosmological solutions jTlj). 

The energy-momentum tensor of (f> is given by 



(4) 



For this energy-momentum tensor one can write the Ein- 
stein equations. 



8ttG 



T, 



,4 1 



(5) 



in the simplified form 



(6) 



Here the gravitational constant is absorbed via the redef- 
inition of the scalar field 



V4^ 



(7) 



It can be shown that the equations ^ are equivalent 
to the Einstein equations for a perfect fluid, which obeys 
the equation of state of stiff matter. A general perfect 
fluid has as energy-momentum tensor given by 



(8) 



where e is the energy density, p is the pressure of the fluid, 
and u'^ is the normalized 4- velocity vector, u'^u^ = 1. A 
stiff fluid has the equation of state e = p and thus the 
speed of sound in the fluid equals the speed of light ■ 
By making to following associations: 

e^p^d.ad-'a, ^ f'" , (9) 

it was shown that a non-rotating stiff fluid is equivalent 
to a massless scalar field ([7]) [13]. 
Inserting ^ into ^ we obtain 

^ -S' -V^^ , 



Rtt 



S-S" 



-v + v" + 



1 + kz 
k 



1 + kz 



k 



v = 



y = 



(10) 



Rzz — — 



Rz 



l + kz 
k 



S' - V'^ + S" -S^ 4cr'2 , 
S - VV = 4(70-' . 



1 + kz 

With some slight manipulation one can re-write the sys- 
tem of equations for a{t, z), S{t, z) and V{t, z) as 



(7 — cr — 



l-f fcz 



a' ^V-V" 



k 



l + kz 



V = Q 



S' = - 



S = -^^^^iV'V + i&a') , (11) 



2k 



V'- 



+ 4(3-^ + 4o-'- 



Note that the scalar field, cr, and the metric ansatz func- 
tion, V{t,z), associated with the x,y coordinates obey 
the same equation. In the next section we present our 
standing wave solution to this system. 



III. STANDING WAVE SOLUTION 

One can show that 

V{t, z) = A Jq + Lozj cos(wi) , 



(T(t, z) = ^ Jo + ^zj sin(wt) , 



(12) 
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solves the two equations from the first hne of (|lip . where 
Jo is the ordinary Besscl function of zeroth order and A 
is the wave amphtude. In addition plugging these func- 
tions into the second line of gives S — 0, i.e. the 
ansatz function S{t, z) — S{z) is time- independent. This 
roughly corresponds to intuition from the electrodynam- 
ics, where it is only the fields and potentials in the trans- 
verse X, y directions which are time dependent. Finally 
taking V{t,z) and <j{t,z) from (fT2)) . inserting them into 
the last equation of pT]) and integrating gives 



■Ml 



LOZ ) J2 ( "T + 



(13) 



Here C is an integration constant and Ji and J2 are or- 
dinary Bessel functions of first and second order respec- 
tively. 

In order to interpret this solution as a domain wall 
embedded in some scalar field which extends to the right 
hand side (i.e. z > 0) we should recover the pure domain 
wall metric as 2 ^ which is the location of the delta 
function energy-momentum tensor source of the domain 
wall. In other words we want V{t, z) — > and S{z) 
as z 0. In the limit z — > we find for V{t, z) 

y(t,0)^A Jo(^)cos(wt) . (14) 

To have V{t,{)) ^ we need to chose uj/k — z„ where 
z„ is the n*'* zero of the Jq Bessel function. If we take 
the properties of the domain wall as fixed (i.e. k or the 
domain wall energy density/tension is fixed) then the re- 
lationship implies that only certain discrete frequencies 
of oscillation are allowed, namely uj = kzn- Note also 
that the scalar field, a, will also vanish at z = under 
the condition u> = kzn- This condition also makes the 
first and last term in brackets of vanish. To can- 
cel the middle term one needs to choose the integration 
constant as 



C 



(15) 



where again z„ = uj/k. 

Examining the energy density associated with the 
scalar field one can show that its energy density is posi- 
tive definite. Using the expression for a from (fT^ in (|3|) 
we find 

Too = 2do(Tdoa - goo {d^adoo + d^ad^cj) = 

■a , ^ j2 ^^^^^^ COS^{Lut)+ (16) 
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i (^— + ujz^ sin^(cjt) 



> . 



However, the stiff fiuid version of the matter source (i.e. 
([5]) plus has an energy density which oscillates be- 
tween positive and negative values. This can be seen by 



inserting a from (fT2l) into the stiff fluid energy density 
e = d^ad^a. Currently many types of matter are consid- 
ered which violate some or all of the energy conditions 
that were formerly imposed on matter sources in general 
relativity (e.g. dark energy [l3|, phantom energy 15]). 
Therefore one might chose not to worry about the non- 
positive definite feature of the stiff fluid. Alternatively 
one could introduce a second fluid which acts as a cos- 
mological constant and alter the stiff fluid to an almost 
stiff fluid as in [l3| where it was shown that this gave 
a positive definite energy density for the fluid version of 
the matter source. 

The standing wave solution given by ([2]) and p2)) can 
be said to have a single polarization. For traveling waves 
one can write a linearized form of the metric ^ as (see 
section 15.2 of [4]): 



ds^ = dt^ - (1 + fxx) dx^ - 2fxydx dy - 
- il + fyy)dy^-dz^ , 



(17) 



where fxx/ 



,^xyF{u) and u — t ~ z] F(u) can 



be any function, although usually it is taken to be some 
sinusoidal function, e.g. F{u) — cos(u) or sin(u). 

There are two different 'linear' polarizations that are 
possible (i) axx = -ayy and axy = 0, or (u) axy = 

o-yx 7^ and Qxx = o-yy — 0. In the limit of large z using 
the standard expansion of Bessel functions of the first 
kind 











UJZ - 
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i] 











Jn(wz) f 

the ansatz function V{z,t) becomes: 

y(z, t) = AJq + LUz'j COs{LUt) 

« A 



(18) 



\ cos (ujz \ cosiujt) . (19) 

V TTLOZ \ 4/ 



Then expanding exp[±y(z, t)] « 1 ± V{z, t) the dx'^,dy'^ 
part of the metric ^ becomes: 



kz [(1 + V)dx^ + (1 - V)dy^ 



(20) 



where V{z^t) is given by (|19p . Ignoring the overall pre- 
factor of kz, and comparing ((20l) with the dx^,dy^ part 
of the metric pT)) . we find that the solution given by ^ 
with has linear polarization, 



^xx ^A/ ^yy ' ^xy 1 

TT 



(21) 



corresponding to case (i) for the metric pT)) . 



IV. DISCUSSION OF PHYSICAL MEANING OF 
SOLUTION 

In this article we have given a simple example of a 
standing gravitational wave with plane symmetry. This 
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solution is not a vacuum solution but requires the pres- 
ence of a domain wall at z — and a massless scalar field, 
cr, occupying the region z > 0. From (|f 2p one finds that 
the metric function, V{t,z), and the scalar field, a(t,z), 
had the same spatial dependence but their time oscilla- 
tions where tt/2 out of phase. One could view this as 
the energy of the oscillation passing back and forth be- 
tween the scalar field and the wave (time-dependent) part 
of gravitational field, as embodied in the metric ansatz 
function V{z,t). Since the field energy of the scalar field 
can be localized this phase connection between the oscil- 
lations of the scalar field and the gravitational field sug- 
gests that for this solution one might be able to define 
a local gravitational field energy - at least for the wave 
(time-dependent) part of the solution. As mentioned in 
the introduction a definition of local gravitational field 
energy is strictly impossible because of the equivalence 
principle. Nevertheless the phase relationship between 
the scalar field and the dx^ , dy^ components of the met- 
ric (exp[±y (z, t)] factors) is suggestive of the wave part of 
the gravitational field carrying a localized energy. Note, 
since this solution is not asymptotically flat and the do- 
main wall is not a localized source, one can not define 
global gravitational field energy via surface integrals over 
effective energy-momentum tensors as discussed in sec- 
tion f4.3 of reference [3] or section 20.5 of reference [1]. 

Absent the matter (i.e. the scalar field) needed to 
support this gravitational standing wave, there is a lot 
in common between the present solution and the simple 
electromagnetic standing wave between two infinite con- 
ducting planes. One plane is obviously the domain wall 
at z = but where is the other plane? A possible answer 
lies in the static, Newtonian potential 



(22) 



^(z) = 2 [500(2:) - .900 (^ = 0)] 



In the limit of large z the ansatz function S{z) from (fT3l 
takes the form 



^ 

^ ^ 27T 


cos^ 


— COS (^ujz 


vr 

~ 4 


Aujz 




TT 





- 2 cos 



Stt 

T 



cos UJZ — 



(23) 



In the first step we have used (|I8p and in the last step we 
have combined all the cosines using standard identities. 
Using these results the Newton potential takes the 
asymptotic form 



$(z) 



exph£(J 
2y/l + kz 



exp [Aloz/tt] 
^fkz 



(24) 



This exponentially increasing, static potential may be 
thought of as trapping the orthogonal oscillations of the 
gravitational field (i.e. those associated with the metric 
components gxx^Qyy and the ansatz function Viz^t)). In 



other words it is the increasing static gravitational field 
associated with ijoo, which traps the oscillatory parts of 
the gravitational field associated with gxxiQyyi thus act- 
ing as the second (soft) plane in conjunction with the 
(hard) plane of the domain wall at z = 0. Also note that 
just as electromagnetic standing waves have discrete fre- 
quencies, so too the plane gravitational standing waves 
have discrete frequencies. For the the present system one 
finds w = fczn, where z„ are the zeros of the Jq Bessel 
function and fc is the energy density/tension of the do- 
main wall. 

The solution presented here also has similarities to the 
neutron trapped in a gravitational field experiments of 
Q . In these experiments cold neutrons were trapped be- 
tween a reflecting surface on the surface of the Earth 
(the domain wall in this paper) and at their upper ex- 
tent by the Earth's gravitational field (the z ^ 00 grav- 
itational field of the domain wall and scalar field in this 
paper) . The energy of the neutrons in these experiments 
was quantized as is the case of the gravitational standing 
waves of the present paper. Note that for neutrons there 
is a well defined concept of local energy. 

One can combine the various asymptotic forms of the 
ansatz functions given in (flQl) and ((23)) to write the metric 
© as: 



ds^ 



exp[j4ajz/7r] 



fcz 



V kz 
1 + A 

\-A 



{dt^ - dz^) - (25) 
cos {ijjz — cos(a;i) 

TTWZ V 4 / 



dx' 



■ cos 



{^z — — ^ cos(ci;t) 



Written in this asymptotic form one sees that the metric 
splits into the usual, time-dependent standing wave form 

(i.e. the i^y^^^cos {ujz — -j) cos(a;t) terms) sitting in 
a static background (i.e. the "^pI^-^^/^I terms). This 
split gives further support to the suggestion that for this 
solution the phase relation between the scalar field and 
ansatz function, y(z,i), may allow one to think about 
a local gravitational field energy for the time-dependent 
oscillating part of the metric - but not the static, back- 
ground part of the metric. 

Although this system is not realistic in the sense that 
domains walls are ruled out (in the present era) by obser- 
vational evidence [l^ it is nevertheless an example of an 
exact, plane standing wave solution to the full non-linear 
Einstein field equations, which has many similarities to 
the simple electro-magnetic plane standing wave. Fur- 
thermore there are two features of this solution which 
may point to some kind of localization of gravitational 
field energy: (i) the oscillation with 7r/2 phase difference 
between the scalar field, a, and the transverse compo- 
nents of the metric, gxxiQyy] (ii) the apparent trapping 
of the transverse oscillations of the metric components 
gxx,gyy by the static gravitational field associated with 
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.900- 

One could also use the present 4D solution to study 
similar standing wave solutions in 5D or higher brane 
world models proposed in (isj . Investigations in this di- 
rection arc underway p^ . 
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